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Abstract 



We examine the effect of one loop quantum corrections on the formation of 
nontopological solitons in a strongly coupled scalar-fermionic Yukawa theory. 
The exact one fermion loop contribution is incorporated by using a nonlocal 
method to correct the local derivative expansion approximation (DE) of the 
effective action. As the Yukawa coupling is increased we find that the nonlocal 
corrections play an increasingly important role. The corrections cause the 
scalar field to increase in depth while maintaining its size. This increases 
the energy of the bag configuration, but this is compensated for by more 
tightly bound fermionic states with lower energy. In contrast to the semi- 
classical picture without quantum corrections, the binding energy is small, 
and the total energy scales directly with the Yukawa coupling. This confirms 
the qualitative behavior found in earlier work using the second order DE, 

although the quantitative solutions differ. 
ll.10.Lm, 03.70.+k 



Typeset using REVTgX 



In theories that support static nontopological soliton solutions at the semi-classical level, 
it is of interest to ask what effect quantum corrections may have on the form of the solution. 
In determining these solutions, one must face the issue of self-consistency — solving a set 
of coupled nonlinear differential equations to find the field configuration that minimizes the 
energy functional. This makes it useful to have a local expansion of the effective action 
to model these vacuum effects. One such local method is to expand the one loop effective 
action in momentum space about zero momentum. For nontrivial background fields this is 
an expansion in increasing orders of background field derivatives, the so called derivative 
expansion (DE). Such an expansion has been considered in theories of finite nuclei Jl|, QED 
0, and Friedberg-Lee type soliton models || ||. 

Of considerable interest is whether the DE converges rapidly enough to make it a useful 
tool. The nature of the convergence has been studied by several authors In a recent 

study [|1(J, we have devised a method in 3 + 1 dimensions that enables the DE contribution to 
the field equations to be corrected for cases where it does not converge. This allows the one 
loop fermionic vacuum scalar density to be included exactly. The method relies on noticing 
that a DE of the fermion Green function S(x, x; u) works well for large loop energies uj 
and for higher partial waves. This method involves making a partial wave expansion of the 
scalar vacuum density, and correcting smaller partial wave DE contributions by an exact 
evaluation of the fermionic Green function. Details of the method and a discussion of its 



reliability can be found in [10 



In this article we consider the effect that vacuum corrections can have on the self- 
consistent solutions of a soliton model. For simplicity, we consider the nontopological soliton 
model of Bagger and Naculich These authors solve for bound states of N fermions with 
a Yukawa coupling to a dynamic scalar field, while including the effects of the one loop 
fermion vacuum using the DE to second order in derivatives of the scalar field. However, 
it is not clear a-priori that the DE is applicable for this model. This is because the DE is 
an asymptotic expansion in inverse powers of mR, where m is the fermion mass and R is 
typically the soliton size or surface thickness. Thus, it is only appropriate for sufficiently 



large solitons with R ^> 1/m. Hence the applicability of the DE will depend on whether 
1/m, or some other scale in the Lagrangian, is instrumental in determining the size of the 
soliton self-consistent solution. 

The choice of Lagrangian density is 0] 



TV 
i=l 
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There are N flavors of fermions in this model, and we work in the limit of large N. The 
utility of the large N parameterization is to validate the semi-classical treatment of the 
scalar field. In particular, the scalar loop contributions are suppressed by a factor 1/N. 
The scalar field has a nonzero vacuum expectation value (VEV), <fi v = ((f)) = yNv, and the 
fermion mass is identified as m = gv. Note that due to the presence of a Landau pole ||, 
only models with g < 30 are physically acceptable. 

To examine the model, we follow Bagger and Naculich || and consider the simple case 
of solitons with iV fermions that all appear in the lowest single particle energy state. The 
fermionic wave function then has the form 
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( iG(r) ^ 
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where \ is a Pauli spinor. By rescaling the scalar field to <p = g<p/vN, the iV dependence 
in the Lagrangian density will contribute simply as an overall factor N. The effective action 
to one loop order in fermionic fluctuations is r e fj[yj] = / d 4 x C + r vac [yj], where r vac [</?] = 
— iTr ln(5 ,_1 ((^)5'((^ i; )), and S(f) = — y?) is the fermionic green function operator in 

the background tp. This gives 

Lvactv 5 ] — — i J d A x tr(x| ln(i$ — ip) — ln(i^ — m)\x). (3) 

The ground state configuration for static fields is then the one that minimizes the energy 
functional E tot [ip] = -T cS [ip}/ J dt. 

After some standard manipulations of Eq. (|3]), the vacuum energy can be written as [J5|.|7| 
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E vac I 



= — — / d 3 x / du utT[j S(x,x;iu)], (4) 
where 

S(x, x; iuS) = (x\- , 1 — r^:\ x ) (5) 

is the fermion Green function in the coordinate space representation. The scalar vacuum 
density is obtained from the energy by taking the functional derivative 5E vajC /5^p, giving 

Pva, c (x) = -— [ dujtv[S(x,x;iuj)]. (6) 

One can derive the DE expressions for E vac and p vac by rewriting the Green function as 



S{x,x;iu) = f —^-(x\p)- (p|ac) 

J (27r) 3 t>y uj - 7 ■ p - <p{iV p ) 

d 3 p 1 



(7) 



(27r) 3 i'jou — 7 • p — tp(x + iV p ) ' 

and expanding about (p(x). Of course, the above expressions for E vac and p vac need to be 
regularized and renormalized to make them finite. Dimensional regularization can be used 
with the appropriate renormalization conditions given in ||. 

It is useful to rescale all fields and variables in terms of the fermionic mass scale m, which 
we can set equal to 1. The effective energy of the system can then be written as an energy 
per fermion in units of m — gv: 

E tot [if] = £ sca l ar [<f] + € + E vac [if] , (8) 
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(9) 

+ O(0 4 ). (10) 



The DE expression for the vacuum energy has been expanded to second order in derivatives 
of the scalar field The fermion energy e is found by solving the eigenvalue equation 

(— «7 ■ V + cp)ip = ejoip- Minimizing the energy functional, -EtotM; and using Eq. (|2]), allows 
us to write the field equations as 
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dr r 

i 

-dn- 
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subject to the normalization constraint 



dG G 

= - +e + ^F, 11 
dr r 

,,F ' {e~9)G, (12) 



u 2 

VV + ^(<P - ¥> 3 ) = /(Pval + Pvac), (13) 



The source terms are 



/ dr{G 2 + F 2 ) = 1. (14) 

M 



P^ = ^r 2 {G 2 -F\ (15) 

Pvac=P v L a D c A +P v D a ^ +0(9 4 ), (16) 
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DE,2 = J_ 
/^vac g7r 2 



^ 3 ln(^ 2 )-^ 3 + ^, (17) 

(18) 



Pvac 4?r 2 

-(Vv.) 2 + ln(^ 2 )(VV) 



As written, these equations facilitate trying different levels of approximation to the 
fermionic vacuum in the self-consistent solution. Setting p vac = corresponds to the semi- 
classical approximation (CL), where vacuum corrections are ignored. Keeping the first term 
in ( ]T6| ) will be referred to as the local density approximation (LDA), as this term gives the 
exact one loop result for a spatially uniform scalar field. Including the second term in ( |T6"D 
will be referred to as derivative expansion (DE) approximation. The density p vac has only 
been expanded to second order — the same order as the differential equation ( |l~3l) for (p. In 
our experience, attempting to include terms with higher order derivatives in self-consistent 
calculations tends to make the solutions numerically unstable. We discuss below a method 
for evaluating p vac exactly. 

It was noted in || that for g — > oo, the soliton solution satisfies a constraint that becomes 
independent of g. To see this we note that all source terms enter with the same power of 
g in (pl|-|l~8"|) , and Eqs. ( |TT] ) and ( jT2P are (^-independent, so with the above notation this 
constraint takes the simple form 
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Pvai(r) + p V ac(r) = 0{l/g 2 ). (19) 

Thus a large amount of cancellation is required between these source terms for large g. 
This restricts the form of the solutions, but does not necessarily enforce an asymptotic g- 
independent shape, as claimed in ||. The reason is that a family of different functions may 
exist that satisfy Eq. ([19]), so that the true solution is determined by the 0(1/ g 2 ) parts. 



The field equations ( |TT| - |T3| ) agree with those of Bagger and Naculich || . Although these 
authors checked the convergence of the DE by examining the relative size of the fourth order 
terms in the expansion, they did so only for the energy functional using a fixed scalar profile. 
The convergence, however, may be quite different for terms in the dynamical equations of 
motion. To see this, consider a fixed background field of the form 



v(r) = i - jy + + / e L • (20) 



At values of a = 0.5, b = 0.16, and / = 0.8 we obtain a scalar field that roughly corresponds 
to the quantum soliton solution at g = 25 found by Bagger and Naculich . Table | shows 
the DE contributions to the energy and the density up to fourth order. Expressions for these 
fourth order contributions can be found in the Appendix. It can be seen that although the 
fourth order terms do not make a substantial contribution to the energy, this is not the case 
for the density (which is shown for r = in the table). For the self-consistent solution the 
convergence is worse as g increases and the size of the scalar bag solution shrinks. 

Our purpose therefore, is to reconsider this model and make use of the correction method 



described in Jl0| to account for the fermionic vacuum in an exact manner. This involves 
extending a scheme discussed by Wasson and Koonin |5]|| for calculations in one spatial di- 
mension. In three spatial dimensions, an exact, or "brute-force" calculation involves making 
a partial wave expansion of the Green function ([5]), 

S(x, x'- iu) = — l Y d S K (r, r'- iu) ® y Rm (x)yl m (x'). (21) 

' ^ K,m 

Each partial wave Green function obeys the equation 
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S K (r,r';iu) =6(r-r'), 



(22) 



which can be solved numerically [Q. A similar partial wave expansion can be made of the 
DE approximation to the Green function, Eq. (|7]), giving @. The sum of all partial 
wave DE contributions to the density will reproduce the expressions ( p^ - |T8| ). 

Wasson and Koonin |5],|9| pointed out that the DE approximation to the Green function 
works well for large loop energies uj (and similarly also for large angular momenta k in 
three spatial dimensions). Therefore one can use the DE as a sophisticated extrapolation 
procedure for accelerating the convergence of the brute-force method. The exact Green 
function contribution can be calculated up to some K max , and the DE can be used to calculate 
the contribution for the remaining partial waves. An equivalent method is to use the full DE 
result, and correct the low energy and low partial wave terms using the difference between 
the exact and DE expressions for S K . This improves the convergence of the energy integral 
in (|6|). Hence we put 



Pvac Hvslc 



+ P^ 2 + APc 



(23) 



where the density correction is given by 

oo 

Apcon#] = J2 Jim Ap K (A K ,r), 



where 



Ap K (A K ,r) = - / du 

7V Z JO 



\x tr[S K (r, r- iu)) + T\<p - V (VV) 
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(26) 
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ip = (p(r), and i K , k K are the modified spherical Bessel functions of order 



k. The utility of this method is that one only needs to include in Ap corr as many partial 



waves as are needed to achieve the desired accuracy. Also, Ap corr is a finite quantity that is 
independent of renormalization. The renormalization counterterms only appear in the DE 
expressions ( |T7| - |T8|) . The details of deriving this correction are described in |T0| . 

Numerically, it is useful to be able to treat our equations entirely as a boundary value 
problem. This can be done by treating e as a field, and also introducing a field \ f° r the 
auxiliary equation QUI) , so that 



dr dr 



The boundary conditions for the entire system of equations are then 



d(p(r) 



dr 
F(r) 



rG(r) 



x(o) 



d(p(r) 



dr 



Ai + i ) [1-<p{R)\, 



b(o) 



G(R) 

x(R) 



g(g) - e 
ifi(R) + e 



1/2 



(28) 



(29) 



where R is large compared to the length scale of the problem. To solve our equations we 
make use of the program COLNEW by Ascher and Bader |13 |. Unfortunately, when the 
correction (|25| ) is included we cannot simply use the COLNEW routine as written, since 
our equations are actually integro-differential equations. To facilitate this, COLNEW was 
modified so that the scalar field solution and its derivative can be extracted at intermediate 
stages of the calculation to evaluate the density correction. The calculation is then internally 
iterative, with Ap corr treated as a source term that is reevaluated as needed with each internal 
iteration of the COLNEW code. 

We now examine the results when using different levels of approximation to the vacuum. 
We use the schemes discussed above in the following order: CL, LDA, DE, and EX. Here EX 
will refer to the exact solution, which includes Ap corr (|24]). Our goal is to quantify the effect 
of this correction on the self-consistent solution for various couplings. The choice g = 10, 
fi = 1/10 will serve as an example. (As the energy scale m already includes an implicit 
dependence on g, setting [i = 1/g simply fixes the scalar field mass at the VEV value.) 
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For the CL, LDA , and DE approximations the self-consistent solutions for ip are shown 
in Fig. |I[ We see that adding the LDA term reduces the depth and width of the scalar field, 
so that a zero no longer appears. The total energy of the solution, E tot = .960m (LDA), has 
dramatically changed from that of the semi-classical approximation, E tot = .620m (CL). This 
behavior, which was identified by H at the DE level, occurs even before the derivative terms 
are included. However, for g > 10.4 the LDA approximation gives no solution. Looking at 
the expression for the LDA contribution in (|T7|) , we see that this function is odd with respect 
to the scalar field. Increasing the coupling drives the scalar field deeper, so that near g = 10.4 
the solution becomes negative at r = 0, and the LDA contributes with the opposite sign. 
At this point the equations no longer support a solution as it is impossible to satisfy the 
constraint ([19]). However, when the DE terms are included the vacuum density is not forced 
to change sign as g is increased, and a soliton solution exists. In Fig. [I] we see that the shape 
of the scalar field for g = 10 is even shallower for the DE than for the LDA. An interesting 
feature to note is that under the DE the total energy, E tot = .967m (DE), is relatively 
unchanged from the LDA result. This occurs even though the Fermi level has changed from 
e = .863m (LDA) to e = .908m (DE), and is made possible due to a corresponding decrease 
in the energy contribution from the scalar field, -E sca iar + -^vac- 

Now examine what happens when the density correction terms are included. In Fig. |2| 
the self-consistent scalar field solutions are shown for calculations including an increasing 
number of partial wave corrections. The exact vacuum density favors a deeper scalar field 
than the DE. The fifth partial wave has not been included in the figure since the field in 
this case is found to be indistinguishable from that with four partial waves. Looking at the 
sequence of scalar fields in Figs, p] and ^| we see a manifestation of the fact that terms in 
the derivative expansion display oscillatory convergence ||12|| . Note that after including the 
k = 1 term in the correction fl25|) the correct vacuum density is specified at the origin as 
higher terms are zero there. However, the k > 1 terms may still affect the self-consistent 
solution. In fact, we see in Fig. ^|that the scalar field becomes deeper as subsequent terms in 
the correction are included. The corrected solution gives e = .857m (EX), and E tot = .976m 
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(EX). The energy of the scalar field is twice the DE result, and the Fermi level has dropped 
away from m. However, the total energy remains fairly stable, rising only slightly above the 
DE result. 

To get the EX correction it is sufficient to use partial waves up to k = 4, except at large 
g where k = 5 is needed. In Fig. ^ we see that for other values of g there is also a decrease 
in the depth of the scalar field using the exact calculation. The exact results agree with 
the DE calculation for small coupling as expected. As g is increased the depth of the bag 
decreases. For g > 20 values of the wave function (p(r) near r = are numerically uncertain 
by about 0.02. This uncertainty is a symptom of the constraint (|i~9|), but in no way affects 
the value of the energies. In Fig. ^ we show how including the correction affects the energies 
at different g. ^From the fermion energy, e, we see that the corrections affect the solution for 
all g > 5. As the coupling increases we have found that not only does the width of the bag 
shrink, but the depth increases. Since this leads to larger derivatives, the energy of the scalar 
field increases. Correspondingly, fermions may be more tightly bound in the bag so that 
the Fermi energy level decreases. When the fermionic and scalar field energies are added to 
obtain the total energy, the sum gives an answer that is closer to the DE result than either 
of the components separately. This occurs at all g. In fact, it is surprising to find that the 
corrected total energy remains fairly constant as the coupling is increased, exhibiting nearly 
the same scaling as the DE total energy. Recall that our energy is scaled in terms of the 
fermionic mass m = gv, so that relative to v the energy scales directly proportional to the 
coupling. Thus the claim that such scaling behavior is universal |14| is supported by our 
calculations. 

In summary, we have found a way of correcting the DE approximation for fermionic 
vacuum fluctuations while still demanding a self-consistent solution. There remain several 
ways in which the results of the calculation done here may be made more general. To include 
different fermionic energy levels it would likely be more economical to use the COLNEW 
routine to solve only the scalar equation, while solving the eigenvalue problem for the relevant 
fermionic wave function components using standard Runge-Kutta techniques. If one takes 
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for granted some validity in the one loop approximation for finite N and large coupling, it 
should also be possible to include the scalar fluctuations at the one loop level. It is expected 
that these fluctuations would contribute a source density term with the opposite sign to the 
fermionic fluctuations 0. However, it is likely that a similar density correction calculation 
would be necessary for the scalar loops to ascertain that this contribution to the true scalar 
vacuum was also included correctly. 
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APPENDIX A: FOURTH ORDER DERIVATIVE EXPANSION OF THE 



The expressions for the energy and density at fourth order, which were used in Table |, 

are 



ACKNOWLEDGMENTS 



ENERGY AND DENSITY 




DE,A 




1607T 2 
1 



1 



/ 



X 



(<9V) 2 _ ll(d aV ) 2 (d 2 v ) 11(^) 4 - 

<p 2 9 V 3 i8<^ 4 y 

25(d a d 2 ip)(d a ip) ll(d 2 d aV )(d a ip) lQ(d 2 V ) 2 



(Al) 



vac 



(A2) 



Covariant notation is used for convenience. 
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TABLES 

TABLE I. Contributions from different orders of the DE series for the energy and the density. 
Results are shown for a fixed background field (p(r) with a = 0.5, b = 0.16, and / = 0.8. Units are 
given in terms of the scale m, as indicated. 





LDA 


2nd order 


4th order 


E vac (10~ 2 m) 


5.909 


2.001 


-0.636 


Pvac(r = 0) (10- 3 m 3 ) 


-5.109 


-0.936 


-3.732 
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FIGURES 




FIG. 1. Self-consistent scalar field solutions of the Bagger-Naculich model (g = 10, fi = 1/10) 
for semi-classical (CL), local density (LDA), and derivative expansion (DE) approximations to the 
vacuum densities. 
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FIG. 2. Self-consistent scalar field solutions with increasing number of terms in the partial wave 
correction series, g = 10, n = 1/10. 
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FIG. 3. Depth of the self-consistent scalar field solutions ip(r = 0) for different values of the 
coupling g, and n = 1/g. Results are shown for the semi-classical (CL), derivative expansion (DE), 
and exact (EX) solutions. 
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FIG. 4. Self-consistent fermion (e) and total soliton (E to t) energies for different values of the 
coupling g, with fi = 1/g. 
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